Developing methods for identifying the inflection point of a convex/concave curve 
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We are introducing two methods for revealing the true inflection point of data that contains or not 
error. The starting point is a set of geometrical properties that follow the existence of an inflection 
point p for a smooth function. These properties connect the concept of convexity /concavity before 
and after p respectively with three chords defined properly. Finally a set of experiments is presented 
for the class of sigmoid curves and for the third order polynomials. 
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INTRODUCING GEOMETRICAL METHODS 

We are going to present some new methods 
for identifying the inflection point of any given con- 
vex/concave or concave/convex curve based on the defini- 
tion and the geometry of it. Before starting it is necessary 
to give some preliminary definitions. 

Preliminary Definitions 

Let a function / : [a,b] -» R, f G C^, n > 2 which 
is convex for x 6 [a,p] and concave for 16 [p, b], p is the 
unique inflection point of f in [a,b] and let an arbitrary 
x <E [a, b]. 

Definition .1 Total, left and right chord are the lines 
connecting {(a, /(a)), (b, f(b))} , {(a, /(a)), (x, f(x))} 
and {(x, f(x)),(b, f(b))} with Cartesian equations g(x), 
l(x) and r(x) respectively. 

Definition .2 Distance from total, left and right chord 
are the functions F, Fi,F r : [a, b] — > R with: 



F(x) = f(x) - g{x) 
F l {x) = f(x)-l(x) 
F r (x) = f(x)-r(x) 



the algebraic surfaces: 
si(a,x) 
s r (b, x) 



Fi{t)dt 



F r (t)dt 



(4) 

(5) 
(6) 



Definition .4 x-left (xi) and x-right (x r ) are such that: 

(7) 

(8) 



xi = argmin {si (a, a;)} 
x£ [a,fc+5i] 



with 61,62 > taken as small as necessary for x\ , x r to 
be unique unconstrained extremums in the corresponding 
intervals. 

We call standard partition (SP) the strictly sorted grid 
of points, not necessary equal spaced: 

Xi, i = 0, 1, . . . n, a = xo < xi < . . . < x n = b (9) 

The corresponding (xi, fi) data produced from f by the 
no error process: 



fi = f(xi),i = 0,l,...n 



(10) 



If errors occur then we have the (xi,<f>i) errored data 
produced from f by the process: 



<Pi = f(Xi) + €i 



iid(0, a 2 ) i = 0,1,... n (11) 



Our analysis is focused on uniform distributions (e, ~ 
U(—r,r)) but it is applicable for every distribution with 
zero mean, for example the normal N(0, a 2 ). If the error 
distribution is not a zero mean one, then the results are 
ambiguous. 

Definition .5 For the errored data llll we have that: 



(1) 


$(Xi) = (, 


*)(xi) - g(xi), 


i = 0, 1, . . . ,rt 


(12) 


(2) 


$l{Xi) = 


j>{Xi) - l(Xi), 


i = 0, 1, . . . , n 


(13) 


(3) 


$ r (Xi) = (, 


b(xi) - r(xi), 


i = 0, 1, . . . , n 


(14) 
(15) 



x r = argmax {s r (b,x)} 

[a— 52,b] 



Definition .6 Function f is called symmetrical around 
inflection point or symmetry around inflection point ex- 
ists when: 

f (p + x) - f (p) = / (p) - / (p-x) , Vx G R (16) 

or 

f(p + x)+f(p-x)-2f(p)=0,Vx£R (17) 

Definition .7 Function f is called locally (e, 6) asymp- 
totical symmetric around inflection point or local e — 6 
asymptotic symmetry exists when: 

\f(p + x) + f(p-x)-2f(p)\<e, Vxe(p-6,p + 6) 

(18) 
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Definition .8 For a function f : [a, b] — > R we have that: 

1. Data symmetry w.r.t. inflection point exists if p — 
b = a — p 

2. Data left asymmetry w.r.t. inflection point exists if 
p — b < a — p 

3. Data right asymmetry w.r.t. inflection point exists 
if p — b > a — p 

Definition .9 A function f : [a, b] — > R is called totally 
symmetric or total symmetry exist, if it is symmetric 
around inflection point p and also exist data symmetry 
w.r.t. p. 

Definition .10 For every subsequent Xi < Xj the ele- 
mentary trapezoidal estimation holds: 



fx' f( x ) dx ~ T i,j(f, x i, x j) 
/(x,)+/(x,) , X 
2 



(19) 



And for every standard partition the total trapezoidal 
estimation holds: 



J a f(x)dx « T n+1 (f,a,b) 

zU=0 T i,i+l(f, x i, x i+l) 



(20) 



The Extremum Surface Method 

We can prove that: 

Lemma .1 x-left (xi) and x-right (x r ) are the points 
where left and right chord respectively are tangent to the 
graph G(f). 

Proof 



1. Let xi < p be the first point where l(x) cuts G(f) 
from left to right. The function is convex for 
x £ [a,p], so G(f) is always below the left chord, 
thus giving a negative value for si (a,xi), which is 
increasing in absolute value as xi departures from 
a to the right until point p. After inflection point 
f is increasing, so it is possible to continue adding 
negative values of surface until the point x* £ (p, b] 
for which G(f) and l(x) have one only common 
point. 

If wc continue beyond this point, then 
G{f) and l(x) have again two common 
points (x\, yi), (x 2 , 2/2), x\ < X2 such that 
si(a,xi) < and si(x\,X2) > 0, thus we have 
started adding positive values to the si (a, x) 
and this leads to a raise of the total value 
si (a, x) = si (a, x\) + si (xi,X2)- 
So, function si (a, x) is decreasing for x £ [a, x*] 
and increasing for x £ [x* , b] , thus x* > p is a local 
minimum and we call it xi. 



2. Let x r > p be the first point where r(x) cuts G(f) 
from right to left. The function is concave for 
x £ [p, b], so G(f) is always above the right chord, 
thus giving a positive value for s r (b,x r ), which is 
increasing as x r departures from b to the left until 
inflection point p. After that point f is still above 
the right chord until the point x* £ [a,p) for which 
G(f) and r{x) have only one common point. 
If we continue again beyond this point, then 
G(f) and r(x) have again two common points 
(xi,yi), (£2,2/2), xi > x 2 such that s r (b, x\) > 
and s r (xi, X2) < 0, so we have started adding neg- 
ative values to the s r (b, x) and this leads to a re- 
duction of the total value s r (b,x) = s r (b,Xi) + 

Sr (X1,X 2 ). 

So, function s r (6, x) is increasing for x £ [a, x*] 
and decreasing for x £ [x* , b] , thus x* < p is a local 
minimum and we call it x r . 

Corollary .1 

*«= arg (/( a ) = MlMl (21) 
xela.b+8!] I x — a ) 

x r= arg f/(«) = %M) ( 22 ) 



iE[a- S 2 ,b] 



b- 



with Si, 5 2 > taken as small as necessary for xi , x r to 
be unique unconstrained solutions in the corresponding 
intervals. 

Corollary .2 Let a function f : [a,b] — > R, f £ 
, n > 2 which is convex for x £ [a,p] and concave for 
x £ [p, b] . Then we have one of the following possibilities: 

1. If Xi,x r £ [a, b] then a < x r < x\ < b 

2. If Xi £ [a, b] then xi > b 

3. If x r ^ [a, b] then x r < a 

We define the next theoretical estimator of the inflection 
point: 

Definition .11 The theoretical extremum surface esti- 
mator (TESE) is 



, x h x r £ [a, b]' 
{ $±p* , xi > b 



(23) 



x r < a 



Lemma .2 If the mesh A(n) of the standard partition is 
such that: 

lim nX(n) 2 = 

n— 1 00 

then T n+ \{(j),a,b) is a consistent estimator of the true 
value ofT n+ i(f,a,b). 
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Proof 

For every subsequent X{ < Xj+i the elementary trape- 
zoidal estimation is: 



and the total variance is: 



Ti i+i (0, Xi,Xi + i) — 



{Xi) + 



-0(^+1 ) 

(24) 



(25) 



Taking the expected value we obtain: 

E (Ti,i+i((f>,Xi,x i+1 )) = 
= ^±^ E (</>&)) + ^±^E (<l>(x i+1 )) 
= ^± f{xi) + ^±f {x . +l) 

= r»,t+i(/, Xi,x l+ i) 

so from the linearity of expected value we have also that: 
E (T n+ i(</>, a, b)) = 

Thus our estimator is unbiased. 

We continue by computing the variance of the elementary 
trapezoidal estimation: 



V (Ti t i + i((j), Xi, Xi+i)) 

2 / \ 2 

x i + 1 —Xj ■ 
2 

4 2 4 



V(0(x i ))+(s±pi)V(0(x i+1 )) 



(27) 

We have two cases. 

If standard partition is equal spaced, then Xi+\ —Xi = 
and we obtain: 



V (Ti !i+ i(4>,Xi,x i+1 )) = 



2n 2 



a 2 (28) 



Let' s compute now the variance of estimator 
T n +i(<j>, a, h): 



V 



(T n+ x(4>,a,b)) = V (Y^ = o T i!i+1 ((f>,Xi,x i+ iy\ 
= nV (T iti+1 ((f),Xi,Xi + i)) 



- „ (6 ~ a) V 

= 

2n 



so clearly it holds: 



lim V (T n+1 (4>,a,b)) = lim 



(29) 



2n 



a 1 = 



For the second case, if standard partition is not equal 
spaced then the mesh or norm of the partition is 



A(n) = max (x%+\ — xA 

i— 0,...,n— 1 



Then it is easy to show that: 



V(T iti+1 (ct>,Xi,x i+1 )) < ^-a 2 (30) 



V {T n+1 (<f>,a,b)) < — n\(nf 







(31) 



from our hypothesis. 
So the estimator is consistent. 

Now we are able to compute using our trapezoidal rule 
l.lOl data estimations for si(xo,Xj) and s r (x n ,Xj): 

Definition .12 

si,j+i(xo,Xj) =T j+ x(<f>i,x ,Xj) (32) 

^r,n— j + 1 [Xn , X j ) — T n —j-\-i Xj , X n ) (33) 

(34) 

It is time to define the next data estimators for xi,x r . 
Definition .13 xu Xr ar e such that: 

Xi = x k (35) 
ji = argmin{sij + i (x a ,Xj)} 

Xr = x jr (36) 

j r = argmax{s r . n -]+i (x n ,Xj)} 
je[o,«-i] 

(37) 

We define now the errored data estimator of the inflection 
point: 

Definition .14 The data extremum surface estimator 
(ESE) is 



XS 



Xi + Xr 



(38) 



Lemma .3 The ESE is a consistent estimator of TESE 
with all relevant integrals calculated via trapezoidal rule. 

Proof 

We have proven in Lemma[2]that trapezoidal rule for the 
errored data gives a consistent estimator for the trape- 
zoidal estimation of the actual data, thus xi > Xr ar e con- 
sistent estimators of the true xi , x r respectively, with rel- 
evant integrals trapezoidal calculated. 
If the interval [a, b] is such that both xi,x r £ [a, b] then 
ESE is a consistent estimator of trapezoidal calculated 
xs = ^- 

If xi > b then recalling Proof of Lemma Q] s; (a, x) is a 

decreasing function, so the minimum xi is achieved when 

Xi = b, the rightmost value of [a, b]. 

If x r < a then recalling Proof of Lemma [J] s r (b, x) is 

an increasing function, so the maximum Xr is achieved 

when Xr = a, the leftmost value of [a, b}. 

So, for every possible case, ESE is a consistent estimator 

of the TESE given by integrals calculated via trapezoidal 

rule. 

Because trapezoidal rule is a good estimation for the true 
values of integrals, we conclude that ESE should also be 
a good estimator of the true value of TESE. 
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The Extremum Distance Method 



3. If Xf,2 ^ [d,b] then xp.2 < b 



Definition .15 xF-left (xf.i) an d xF-right (xf,2) are 
such that: 



xf,i = min F(x) 

x£[a— Si,b] 

xf.2 = max F(x) 

x£[a,b+S2] 



(39) 
(40) 
(41) 



with 61,62 > taken as small as necessary for Xf,i,Xf,2 
to be unique unconstrained extremums in the correspond- 
ing intervals. 



We can prove that: 
Lemma .4 



x F \.2 = arg <j / (x) 

,tG [a— 5± ,6+52] 



f(b) - f(a) 



(42) 



with 6i,6- 2 > taken as small as necessary for Xpi, Xp2 
to be unique unconstrained extremums in the correspond- 
ing intervals. 

Proof 

We have extended the interval [a, b] such that there 
exist (both unconstrained) a local minimum and a lo- 
cal maximum inside. For our convex/concave case let 
p G [a — Si, b + 62], p £ {a, b} is the internal root of 
F(x). Then we have that F(x) < 0, x G [a — Si, p] and 
F(x) > 0, x G [p, b + S2], because function is convex near 
a and concave near b. Thus the local minimum exists at 
[a — Si, p] and the local maximum lies in [p, b + 62] ■ If we 
take the first derivative we have that: 

F\x) = f\x)-g\x) = f'(x)-X = 0^ f\x) = A (43) 



where A = jy "' b _ J n ya> is the slope of the total chord. 



f(b)-f(a) 
b—a 

But the above equation must hold for both local min- 
imum/maximum xpi 2, so it is necessary to hold: 



= /(*„) = M^M (44) 



b — a 

We can also check the second derivative which is 
F"(x) = f"(x) 



(45) 



so it holds F (xpi) = f (xfi) > and F (xp 2 ) = 
/ { X F2) < 0, i.e. we have the correct signs for local 
minimum and maximum respectively. 

Corollary .3 Let a function f : [a,b] —> R, f G 
CW, n>2 which is convex for x G [a,p] and concave for 
x G \p, b] . Then we have one of the following possibilities: 

1. If xf.i, xf,2 G [a,b] then a < xf.i < xf,2 < b 

2. If xf.i ^ [1, b] then xp.i < a 



We define the next theoretical estimator of the inflection 
point: 

Definition .16 The theoretical extremum distance from 
total chord estimator (TEDE) is such that 



xd 



xfa + X F ,2 



(46) 



Now we define the data estimators of Xf,i,xf,2- 
Definition .17 



Xf,i = x n 
ji = argmin{<fr(xj)} 

j£[0,n] 



XF.2 = X h 

32 = argmax{<f>(xj)} 

j€[0,n] 



(47) 
(48) 



The data extremum distance from total chord estimator 
(EDE) is 

Definition .18 

Xf.i + Xf.2 ■„ ^ , An \ 
Xd = g l ® Xf - 2 - Xf ' x ' 9 > 

Lemma .5 The EDE is an unbiased estimator of TEDE. 
Proof 

For all 3>(xj),j = 0, 1, . . . , n it holds: 



E($(x j ))=F(x j ) 



(50) 



so if we take the errorcd data instead of the true data it 
has to be also that: 



E min {$(x,)} = min {F(x~)} (51) 



E ( max |$(x,)i = max \F(xj)} 
\je[o,n] 1 J "J j€[0M 



Iterative application of geometrical based methods 

Another very important oppurtunity is the possibil- 
ity of iterations like the well known bisection method in 
root finding. Recall that for a continuous function if 
f{a) f{(3) < then exist £ G (a, f3) such that /(£) = 0. 
Our ESE method always gives an interval that contains 
the true inflection point p or a point close to the edge a 
or b, if data is just convex (or just concave) and inflection 
point does not exist. EDE method also gives an interval, 
although it is more sensitive to errors, so it does not al- 
ways give a point close to a or b, if simple convexity or 
concavity exist. 
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1. ESE iterative method or Bisection-ESE or BESE 

We apply to initial data (xi, <fii), i = 0, . . . , n the 
ESE method and have the th output for ESE 
method: 



TABLE I: Fisher-Pry sigmoid, p=5, total symmetry, n=500, 
no-error 



jr = 170 


jl = 332 


X r = 4.028 


Xi = 5.972 


Xs = 5.0 


j Fl = 155 


j F2 = 347 


Xfa = 3.884 


Xf,2 = 6.1520 


Xd = 5.0 



r,-(0) ,(0)1 (0) _ (0) _ (0) _ x< Q) +x, <0> 

[Jr , J I J, Xr — Xj{0) , Xl ~ x j(°) > Xs ~ 2 

If and only if jj ' > , then we apply again ESE 
for data: 



{Xi,<pi),t =ji 0) ,-.-,Ji 

I si 



(0) 



and obtain the 1 output for ESE method: 



XjM , Xi 



(i) 



(i) 

x J w , Xs 



x^'+x,' 1 ' 



Symmetric sigmoid curves 

We find the points x\ , xgg that give the first 1% and the 
99% of sigmoid' s capacity L because those points have 
economic sense. So our interval [a, 6] is always relevant 
to the interval [xi,xgg]. 



The Fisher-Pry sigmoid curve with total symmetry 
Let' s take the function: 



We continue until jj < jr"' or until 
X^ ~ X^s < e ' w hh e = 10~ 8 to be a good 
tolerance for all examined data. 



,(fe) 



2. EDE iterative method or Bisection- ED E or BEDE 

We apply to initial data (xi,<pi),i = 0,...,n all 
four methods and have the th output for EDE (iff 
xp.2 > x F,i) and ESE methods: 



r.(0) .(Oh (0) _ 

[Jl ,J2 J; Xf,i - X 3 [°> i X F ,2 - x j(°> i X D 



(0) 



Jr ,Jl 



v (0) _ T v (0) _ „ JO) 
i Xr — x j (0) , Xi — X jW , Xs 



x^. 0) +x, (0) 



If and only if j% > then we apply again all 
four methods for data: 



(xi,<f>i),i =j[°\- 



AO) 
1J2 



and obtain the I s * output for EDE (iff Xf\ > Xf\) 
and ESE methods: 



r .(1) (1) _ (1) _ (1) _ 
[Jl iJ2 J: Xf,1 ~ j^ r> ' XF,2 - ^(D: - 2 



,(1) „•(!) 



(1) (1) 
Xr = X (i) , ; = X (l) , x s 



(i) _ x^+x, (1) 



/ (x) = 5 + 5 tanh (x - 5) (52) 

after which has p = 5, L = 10, Xi = 2.7024, xgg = 
7.2976 and examine it at the interval [2,8] in order to 
have data symmetry w.r.t. inflection point. The function 
is also symmetrical around inflection point, i.e. we have 
total symmetry. 

From Corollary [J] we compute x; = 5.970315941, x r = 
4.029684059, x Fl = 3.850750196, x F2 = 6.149249804, all 
inside [2,8], thus all methods are theoretically applicable. 
We first take n = 500 sub-intervals equal spaced without 
error just for checking our estimators. The results are 
presented at Table |TJ 

We observe that X i = 5.9720, Xr = 4.0280, xfi = 
3.8480, XF2 — 6.1520 are very close to the theoretical 
expected values, so we are on the results of Lemma [751 
The absolutely accuracy from the first apply of all meth- 
ods confirms our theoretical analysis. All important lines, 
curves and points are presented at Figure [T] 



f(x)=5+5 tanh(x-5|, inflection point p = 5 


its-5-, 






f 















EXPERIMENTS AND RESULTS 



We design small experiments by taking a suitable 
smooth function of known inflection point p, an inter- 
val [a, 6] that covers all the possible cases p £ [a 7 b} 7 p < 
a, p > b and we add a uniform error £j ~ U{—r,r) via 
the process ITT1 



FIG. 1: (color online, scaled) Fisher-Pry sigmoid with total 
symmetry and without error 



We next add the error term q ~ t/(— 0.05, 0.05) via 
the process [TT] and run our algorithms again. The results 
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TABLE II: Fisher-Pry sigmoid, p=5, data symmetry, [a, b] — 
[2,8], n=500, error r=0.05 



TABLE III: Fisher-Pry sigmoid, p=5, data left symmetry, 
[o, 6] = [4.2, 8], n=500, no error 



jr = 163 


ji = 339 


Xr = 3.944 


Xi = 6.056 


Xs = 5.000 


j Fl = 155 


jp 2 = 349 


Xf,i = 3.848 


X f,2 = 6.176 


Xd = 5.012 



jr = 2 


ji = 156 


Xr = 4.2076 


X! = 5.3780 


Xs = 4.7928 


J>i = 1 


Jf 2 = 234 


Xf,i = 4.2 


Xf,2 = 5.9708 


Xd = 5.0854 



are presented at Table [TTJ 

Again the estimations are close to the theoretically 
expected and both methods gave the true answer from 
the first apply. We present the ESE and EDE inter- 
vals and estimators together with the true function and 
the errored data at Figure [5J where we present all im- 
portant points xi,x r ,XF,i,XF,2 with the relevant tan- 
gent lines and the size of the minimum/maximum of F. 
Due to the total symmetry the (dashed) line connect- 
ing (xf i, f{xF i)), {xf 2, f(xF 2)) passes from the point 
(p,/(p)j- 



TABLE IV: ESE iterations for Fisher-Pry sigmoid, p=5, data 
left asymmetry, [a, b] — [4.2,8], n=500, no-error 



Xr 


Xi 


Xs 


4.8156 


5.3704 


5.0930 


4.8232 


5.0892 


4.9562 


4.9524 


5.0892 


5.0208 


4.9600 


5.0208 


4.9904 


4.9904 


5.0208 


5.0056 



TABLE V: EDE iterations for Fisher-Pry sigmoid, p=5, data 
left asymmetry, [a, b] = [4.2,8], n=500, no-error 



f(x)=5+5 tanh[X-5), inflection point p = 5 











t 


-Id" 501 



-2 2 4 6 



XF\ 


XF 2 


XD 


Xs 


4.5192 


5.4844 


5.0018 


4.9638 


4.7244 


5.2716 


4.9980 


4.9980 


4.8460 


5.1576 


5.0018 


5.0018 


4.9068 


5.0892 


4.9980 


4.9980 


4.9448 


5.0512 


4.9980 


4.9980 



TABLE VI: Fisher-Pry sigmoid, p=5, data left symmetry, 
[a, b] = [4.2,8], n=500, error r=0.05 



FIG. 2: (color online, scaled) Fisher-Pry sigmoid with total 
symmetry and error r— 0.05 



jr = 3 


ji = 149 


X r = 4.2304 


Xi = 5.3248 


Xs = 4.7700 


JFl = 3 


jp 2 = 231 


Xf,i = 4.2152 


Xf,2 = 5.9480 


Xd = 5.0816 



The Fisher-Pry sigmoid curve with data left asymmetry 

We continue with the same sigmoid function, but now 
we form proper our [a, b] to show data asymmetry w.r.t. 
inflection point. 

Let' s take for example [4.2,8]. If we do our theo- 
retical computations we find xi = 5.974322740, x r = 
4.029684059, x F1 = 4.025677260, x F2 = 5.974322740. 
We have that x r < a, so Xr has to estimate a = 4.2 and 
Xs must be close to 4.703504993. Additionally, xp\ < a, 
so xf\ must be also an estimation of a, thus \d must 
lie near the value 5.087161370. It' s time to see if our 
theoretical predictions will be confirmed by experiment. 
We use for comparability the same Standard Partition as 
before and have the output presented at Table IIIII We 
have confirmation of our theory 

It is time now to try iterations based on ESE and EDE 



intervals that contain inflection point and to observe re- 
markable convergence to the real value of p = 5 for both 
methods. We present ESE iterations at Table IIVI and 
EDE iterations at Tabk|V] 

Let' s add the same error term <~ U{— 0.05, 0.05) 
and run our algorithms. The results are at Table IVII 
and clearly we are close enough to our theoretical 
expectations. 

We observe that ESE method did not estimae the in- 
flection point with acceptable accuracy, so it is time to 
run the ESE and EDE iterations. The results, Table IW] 
and Table IVIII1 show a clear improvement of both esti- 
mations. 

All the points of interesting are presented at Figure El 
where we see that interval does not contain both xi , x r 
and xf,i, Xf,2- 
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TABLE VII: ESE iterations for Fisher-Pry sigmoid, p=5, data 
left asymmetry, [a, b] = [4.2,8], n=500, error r=0.05 



Xr 


Xi 


Xs 


4.8156 


5.3248 


5.0702 


4.9144 


5.1576 


5.0360 



TABLE VIII: EDE iterations for Fisher-Pry sigmoid, p=5, 
data left asymmetry, [a,b] = [4.2,8], n=500, error r=0.05 





XF 2 


XD 


Xs 


4.5268 


5.5148 


5.0208 


4.9258 


4.7244 


5.2412 


4.9828 


4.9828 



Non symmetric sigmoid curves 

We continue our study with non symmetric sigmoid 
curves that appear in Economics and other disciplines. 

The Gompertz sigmoid curve 

Let' s examine the function: 

f{x) = 10e- c5e " x (53) 

after in the interval [3.5, 8]. The basic properties are 
presented at IIX1 

It is easy to prove that f is (0.224, 1.0)-asymptotical sym- 
metrical around inflection point, so we can hundlc it sim- 
ilar to a symmetric sigmoid only for a distance of ±1 from 
p = 5. 

We use, for comparison reasons, the same SP with 
500 sub-intervals without error and obtain the Table |X] 
which is absolutely compatible with theoretical predic- 
tions. The ESE iterations are showed at Table Ixl . while 
EDE iterations can be found at Table IXII Wc observe 



f{x)=5+5 tanh(x-5|. inflection point p = 5 




2 4 6 8 10 12 



FIG. 3: (color online, scaled) Fisher-Pry sigmoid with data 
left asymmetry and error r=0.05 



TABLE IX: Basic properties of a Gompertz sigmoid 



L — 10 


xi = 3.472820374 


x g9 = 9.600149227 






xs 


4.138928270 


5.887451706 


5.013189988 


Xf,1 


Xf,2 


XD 


4.095750735 


6.290768183 


5.193259460 



TABLE X: Gompertz sigmoid, p=5, asymmetry, n=500, no- 
error 



jr = 72 


3i = 266 


Xr = 4.139 


Xi = 5.885 


Xs = 5.012 


JFi = 67 


JF 2 =311 


Xf,i = 4.094 


Xf,2 = 6.290 


XD = 5.192 



TABLE XI: ESE iterations for Gompertz sigmoid, p=5, asym- 
metry, [a,b] — [3.5,8], n=500, no error 



Xr 


Xi 


Xs 


4.7780 


5.2010 


4.9895 


4.9040 


5.1200 


5.0120 


4.9400 


5.0480 


4.9940 


4.9760 


5.0300 


5.0030 


4.9850 


5.0120 


4.9985 



TABLE XII: EDE iterations for Gompertz sigmoid, p=5, 
asymmetry, [a, b] = [3.5,8], n=500, no error 





Xf 2 


XD 


Xs 


4.4540 


5.6690 


5.0615 


5.0120 


4.6700 


5.3630 


5.0165 


5.0030 


4.8050 


5.2100 


5.0075 


5.0030 


4.8860 


5.1200 


5.0030 


4.9985 


4.9310 


5.0660 


4.9985 


4.9985 


4.9580 


5.0390 


4.9985 


5.0030 



convergence to the real p for all methods used in this 
Chapter. 

We continue with our familiar SP by adding uniform 
error distributed by U(— 0.05, 0.05). Results at Table 
IXTETl while the ESE iterations arc shown at Table IXiVl 
and the EDE iterations at Table IXV1 Convergence to 
the true value occur. All points of interest and data 
arc presented at Figure |H There exists satisfactory 
estimation of the p = 5, so there is no need for further 
investigation. 



Non sigmoid curves 

Our analysis is applicable also to non sigmoid curves, 
not necessary symmetric or with data symmetry. Wc 
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TABLE XIII: Gompertz sigmoid, p=5, asymmetry, [a, b] = 
[3.5,8], n=500, error r=0.05 



jr = 74 


3i = 274 


Xr = 4.1570 


Xi = 5.9570 


X s = 5.0570 


j Fl = 66 


j F . 2 = 319 


Xf,i = 4.0850 


Xf,2 = 6.3620 


Xd = 5.2235 



TABLE XIV: ESE iterations for Gompertz sigmoid, p=5, 
asymmetry, [a, 6] = [3.5,8], n=500, error r=0.05 



Xr 


Xi 


Xs 


4.6340 


5.5340 


5.0840 


4.8590 


5.1560 


5.0075 



TABLE XV: EDE iterations for Gompertz sigmoid, p=5, 
asymmetry, [a, b] = [3.5,8], n=500, error r=0.05 





XF 2 


XD 


Xs 


4.472 


5.642 


5.057 


5.003 



shall proceed making two experiments for a symmetric 
polynomial of r rd order. 

A symmetric third order polynomial with total symmetry 
Let the polynomial function: 

f{x) = -\x* + \x 2 -±x+\ (54) 

We study it at [—2, 7], it has inflection point at p = 2.5 
and we have total symmetry. The interesting points are 
presented here: 









0.25 


4.75 


2.50 


XF,1 


XF,2 


XD 


-0.09807621078 


5.098076211 


2.50 



The SP with 500 sub-intervals without error gives the 
Table IXVll which is absolutely compatible with theoreti- 
cal predictions. There is no need for any kind of iteration, 
because all methods agree. 



f(x)=10e ee . inflection point p = 5 




2 4 6 8 10 12 



FIG. 4: (color online, scaled) Gompertz asymmetrical sigmoid 
with error r=0.05 



TABLE XVI: 3 rd order polynomial, total symmetry, p=2.5, 
n=500, no-error 



jr = 126 


3i = 376 


Xr = 0.25 


Xi = 4.75 


Xs = 2.50 


j Fl = 107 


j> 2 = 395 


Xf.i = -0.092 


X f,2 = 5.092 


Xd = 2.50 



TABLE XVII: Symmetric 3 rd order polynomial, p=2.5, 
n=500, error r=2.0 



jr = 115 


ji = 375 


Xr = 0.052 


Xi = 4.732 


Xs = 2.392 


j Fl = 105 


j F2 = 375 


Xf,i = -0.128 


X f,2 = 4.732 


Xd = 2.302 



TABLE XVIII: ESE iterations for 3 rd order polynomial, p=5, 
total symmetry, n=500, error r=2.0 



Xr 


Xi 


Xs 


1.222 


3.688 


2.455 


1.564 


3.382 


2.473 



The same SP with uniform error distributed by 
U(-2,2) gives the results of Table IXVlTl two ESE itera- 
tions at Table IXVIIII All points and data are presented 
at Figure El 

A symmetric third order polynomial with data right 
asymmetry 

For the same symmetric 3 rd order polynomial [54] we 
change the interval to [—2,8], thus we have data right 
asymmetry now. Our critical points are written here: 

Xr X[ XS 

-0.25 ~ 4.75 ~2^5 

XFS XF,2 Xp 

-0.429732639 5.429732639 2.50 
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t(x)=-1/3 x 3 + 5/2 x E - 4x + 1/2, inflection point p = 2.5 
25 1 ' ' ' T 




FIG. 5: (color online, unsealed) Symmetric 3 rd order polyno- 
mial with error r=2.0 



TABLE XXII: Symmetric 3 rd order polynomial, p=2.5, 
n=500, [-2,8], error r=2.0 



jr = 88 


ji = 338 


Xr = -0.26 


Xi = 4.74 


Xs = 2.24 


M = 88 


jr 2 = 384 


X f,i = -0.26 


Xf,2 = 5.66 


Xd = 2.70 



TABLE XXIII: ESE iterations for symmetric 3 rd order poly- 
nomial, p=2.5, n=500, [-2,8], error r=2.0 



Xr 


Xi 


Xs 


1.46 


3.84 


2.65 



TABLE XIX: Symmetric 3 rd order polynomial, p=2.5, 
n=500, [-2,8], no-error 



jr = 88 


ji = 338 


Xr = -0.26 


X! = 4.74 


XS = 2.24 


j Fl = 80 


Jf 2 = 372 


Xf.i = -0.42 


Xf,2 = 5.42 


XD = 2.50 



TABLE XXIV: EDE iterations for symmetric 3 rd order poly- 
nomial, p=2.5, n=500, [-2,8], error r=2.0 





XF 2 


XD 


Xs 


0.86 


3.84 


2.35 


2.350 



TABLE XX: ESE iterations for symmetric 3 rd order polyno- 
mial, p=2.5, n=500, [—2,8], no-error 



Xr 


Xi 


Xs 


1.3800 


3.8800 


2.6300 


1.8200 


3.0600 


2.4400 


2.2200 


2.8400 


2.5300 


2.3200 


2.6400 


2.4800 


2.4200 


2.5800 


2.5000 


2.4600 


2.5400 


2.5000 



TABLE XXI: EDE iterations for symmetric 3 rd order poly- 
nomial, p=2.5, n=500, [—2,8], no-error 



XFt 


XF 2 


XD 


Xs 


0.82000 


4.1800 


2.5000 


2.5000 



There exist a problem here. Although we have a 
symmetric polynomial, the TESE is not equal to the 
true inflection point. A remedy for this problem for the 
class of 3 rd order polynomials is given with the next 
Lemma. 

Lemma .6 3 rd order polynomial ESE correction. 

Let a 3 rd order polynomial f (x) = ax 3 + Bx 2 +jx + S 
and let p its inflection point. Then it holds exactly that: 

p= \ Xl + \ Xr + \ a + \ b ( 55 ) 

and 

„ 1 1 11, 
3 3 6 3 

is a consistent estimator of trapezoidal estimated p. 



The case of SP with 500 sub-intervals and no error 
gives the Table IXIXl while ESE and EDE iterations are 
presented at Table IXXl and Table IXXll respectively. First 
results are absolutely compatible with theoretical predic- 
tions. For example we are waiting that xs — 2.25 and 
we found 2.24. 

We add uniform error distributed by U{— 2, 2) and we 
have the results of Table IXXII1 one ESE iteration at Ta- 
ble IXXIIII and one EDE iteration at Table IXXIVI All 
points and data are presented at Figure [6] 



Proof 

The inflection point because a ^ is found from the root 

J_ 

3a- 



of the second derivative, i.e. 6 ctp + 2/3 = or p — 



Due to Corollary Q] we have for the Xi that: 

„ ax 3 + Bx 2 +jx — aa 3 — Ba 2 — ^ya 

3 ax 2 + 2/3x + ~f = ■ — 

x — a 

or 

{x - a) 2 (a a + /3 + 2 ax) = 
so the internal solution xi is: 

aa + 8 



xi 



2a 



(56) 
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f(x)=-1/3 x + 5/2 x - 4x + 1/2, inflection point p = 2.5 




FIG. 6: (color online, unsealed) Symmetric 3 rd order polyno- 
mial with data right asymmetry 



For the x r we have similar that: 

a£> 3 + /3& 2 +7& — ax 3 — fix 2 — jx 



3ax z +2f3x + j 



b — x 



or 



{b-x) {ab + j3 + 2ax) = 
so the internal solution x r is: 

ab + (3 

X r — ~ 

2 a 

By adding xi and x r we obtain: 

xi+x r = - -a -b 

2 a 2 

and if we remember that p = — we obtain: 
xi + x r = --a + 3p - - b 



(57) 



or finally 



1 1 1 1 , 

3 3 6 3 



Since we have proven that \i ? Xr are consistent esti- 
mators of trapezoidal calculated values of xi , x r we can 



take a consistent estimation for trapezoidal calculated 
p by replacing the unknown xi , x r with the estimators 

XhXv 

As an example, we come back to the case of 3 rd order 
symmetric polynomial with data right asymmetry. We 
have that a = — 2, b = 8 and from Table IXXIII is \ r = 
—0.26, xi = 4.74, so we have that: 



1 1 

P = g Xi + o Xr 



-a+-b= 2.493333333 
6 3 



which is much closer to the true value of 2.5. 

The above analysis can be extended to every function, 
if we can find analytically a relation between inflection 
point and xi,x r , a, b. 

CONCLUSION 

Starting from the problem of identifying the inflection 
point for the data (ajj, <f>i) i = 0, 1, . . . , n we have created 
two geometrical methods in order to solve this problem: 
the Extremum Surface Estimator (ESE method) and 
the Extremum Distance Estimator (EDE method). 

The methods can be applied iteratively just like bisec- 
tion method and converge to the true i.p. either from the 
first or after a few iterations only. The implementation 
of the methods have been done with Fortran, Maple and 
Matlab. It is interesting to mention that for the prob- 
lem of estimating the inflection point of n = 10000 data 
pairs we needed less than 7 sec CPU time in Fortran 
GNU compiler and a typical Intel Core i5 CPU with 4 
GB RAM memory. 
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